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[57] ABSTRACT 

A real-time contour detector and data acquisition sys- 
tem for an angiographic apparatus having a video scan- 
ner for converting an X-ray image of a structure charac- 
terized by a change in brightness level compared with 


its surrounding into video format and displaying the 
X-ray image in recurring video fields. 

The real-time contour detector and data acquisition 
system includes track and hold circuits adapted to be 
coupled to the video scanner for sampling the bright- 
ness level of predetermined sample points of the current 
video field in the region of the last detected border 
point; a reference level analog computer circuit con- 
nected to the track and hold circuits and responsive to 
the sampling of the brightness level of the predeter- 
mined sample points for generating a reference bright- 
ness level at which the next border point can be de- 
tected; and analog compartor connected to the refer- 
ence level circuit for comparing the brightness level 
along a video scan line in the current video field with 
that of the reference brightness level to detect the next 
border point; a digital processor connected to the com- 
parator and the track and hold circuits and responsive 
to the detection of a border point for enabling the com- 
parator only within an area of the current video field 
covered by a narrow expectation window centered 
about the expected next border position and for deter- 
mining the sample points to be sampled by the track and 
hold circuits; a field memory connected to and accessi- 
ble by the processor for storing the coordinates of the 
detected border points in the preceding two video fields 
to partially define the center of the expectation win- 
dow; and a computer interface connected to the com- 
parators for storing the detected border points in the 
computer. 


8 Claims, 18 Drawing Figures 
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CONTOUR DETECTOR AND DATA ACQUISITION 
SYSTEM FOR THE LEFT VENTRICULAR 
OUTLINE 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 10 
(72 Stat. 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates generally to automatic 15 
detection of structures characterized by a change in 
brightness level compared with its surrounding when 
the information is carried by a video signal, and more 
particularly to the automatic detection of the left ven- 
tricular contour. 20 

2. Description of the Prior Art 

Dynamic measurements of the size and shape of the 
left ventricular cavity and their correlation with simul- 
taneously occurring pressure and flow events are of 
great importance in cardiovascular research. Many 25 
methods are known for measuring cardiac chamber 
dimensions and volume in animals and man. However, 
angiocardiography has proven with time to be the most 
readily available and reliable method for these pur- 
poses. By injecting a suitable radiopaque liquid (usually 30 
Renovist) into the left ventricle, a silhouette of the left 
ventricle can be generated at the face of an image inten- 
sifier with X-rays from an X-ray tube underneath the 
patient. When using the negative of conventional film, 
the X-ray picture shows the left ventricle as a light 35 
structure against a dark background. In the border re- 
gion between the ventricular silhouette and its back- 
ground the so-called ventricular contour can be visual- 
ized. Certain mathematical methods enable the left ven- 
tricular volume to be estimated from this contour data. 40 

The ability to obtain radiographic information con- 
cerning chamber size and shape for clinical or investiga- 
tive use within reasonable periods of time is limited. 
Presently, the acquisition of such data usually requires 
the manual tracing of heart outlines to properly define 45 
chamber borders, even with films of excellent quality. 
The use of more advanced or partially automated sys- 
tems also requires an experienced investigator to define 
ambiguous areas where the diaphragm or ribs cross the 
heart outline or the aorta passes behind the heart image. 50 
Reproducibility often is poor and the task is tedious 
because of the great number of pictures to be processed. 
Since the introduction of the digital computer into the 
cardiovascular laboratory, many aspects of the data 
processing have been automated. However, the crucial 55 
part in the data acquisition and processing is still the 
definition and subsequent storage of the chamber bor- 
ders into the digital computer. 

In the last decade, many attempts have been made to 
automate this procedure and to come up with an auto- 60 
mated border recognizer. Most methods are based upon 
computer processing of the digitized left ventricular 
image, making online and real-time operation impossi- 
ble. 

An electronic system for automatically detecting 65 
structures characterized by a change in brightness level 
when the information is carried by a video signal and 
processing the picture information has been discussed in 
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a paper by Leo P. de Jong and Cornelis J. Slager, “Au- 
tomatic Detection of the Left Ventricular Outline in 
Angiographs Using Television Signal Processing Tech- 
niques”, IEEE Transactions on Biomedical Engineer- 
ing, Vol. BME-22, No. 3, May 1975, pp. 230-237. Start- 
ing from the available angiographic data regarding the 
location and brightness level near an already known 
contour point of the left ventricle for a given video 
scanning line, a prediction is made of the most probable 
location and brightness level of the contour point at the 
next line. The prediction of location, based on the co- 
herence that may be expected between contour points 
on successive lines, is realized by defining a so-called 
expectation window for each subsequent contour point 
and allowing detection only in this expectation window. 
Brightness information obtained at sample points near 
the detected contour is used to predict the reference 
brightness level at which the next contour point may be 
detected. This system is principally based on analog 
circuit design with the expectation and sample-and-hold 
windows determined by time delays generated by 
monostable multivibrators which does not permit a 
relatively easy and accurate correlation of video lines 
and fields. Its other disadvantages include the follow- 
ing. The center of the expectation window on a given 
line is defined as having the same horizontal position as 
the last detected border point on the preceding line. The 
basic limitations of this principle are that only limited 
excursions from the direction perpendicular to the 
video scan lines are allowed and the center of the expec- 
tation window is generally not a good approximation of 
the next border point. The position of the sample point, 
which is determined shortly before the expected new 
border position, is also defined with respect to this cen- 
ter of the expectation window. This means that the local 
direction and movement of the left ventricular outline is 
not taken into account in the positioning of this sample 
point, resulting in a less accurate determination of the 
desired video sample and of the calculated reference 
level. Only two brightness information sample points 
are defined for the calculation of the reference level for 
the next border point to be detected. This does not 
allow for the determination of the sign of the back- 
ground slope in the video signal near the actual border; 
for negative values too high a reference level will be 
calculated according to the formula for the reference 
level. 

BRIEF SUMMARY OF THE INVENTION 

It is therefore one object of the present invention to 
provide an improved contour detector and data acquisi- 
tion system. 

It is another object to provide such a system based on 
a digital approach allowing a relatively easy and very 
accurate correlation of video lines and fields. 

It is a further object to provide such a system having 
potential for online, real-time acquisition of the left 
ventricular shape during angiocardiography. 

The objects of the present invention are achieved by 
a real-time contour detector and data acquisition system 
for an angiographic apparatus having video scan means 
for converting X-ray images of a structure character- 
ized by a change in brightness level compared with its 
surrounding into video format and displaying the X-ray 
images in recurring video fields. The real-time contour 
detector and data acquisition system includes sampling 
means adapted to be coupled to the video scan means 
for sampling the brightness level of predetermined sam- 
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pie points of the current video field in the region of the 
last detected border point; reference level means con- 
nected to the sampling means and responsive to the 
sampling of the brightness level of the predetermined 
sample points for generating a reference brightness 5 
level at which the next border point can be detected; 
comparator means connected to the reference level 
means for comparing the brightness level along a video 
scan line in the current video field with that of the 
reference brightness level to detect the next border 10 
point; processor means connected to the comparator 
means and the sampling means responsive to the detec- 
tion of a border point for enabling the comparator 
means only within a narrow expectation window on the 
next video scan line centered about the expected next 15 
border point and for determining the sample points to 
be sampled by the sampling means; field memory means 
connected to and accessible by the processor means for 
storing the coordinates of the detected border points in 
the preceding two video fields to partially define the 20 
center of the expectation window; and data acquisition 
means connected to the comparator means for acquiring 
the detected border points. 

The foregoing as well as other objects, features, and ^ 
advantages of the present invention will become more L 
apparent from the following detailed description taken 
in conjunction with the appended drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of the contour detector and 30 
data acquisition system of the present invention. 

FIG. 2 illustrates the contour of a left ventricular 
angiogram as determined with the system of the present 
invention. 35 

FIGS. 3 A to 3C illustrate definition of the border 
points with respect to the expectation window. 

FIG. 4 is a diagram for explaining the manner in 
which the center of the expectation window is deter- 
mined by the first order line extrapolation method. 40 

FIG. 5 is a diagram for explaining the manner in 
which the center of the expectation window is deter- 
mined by the zeroth order line extrapolation method. 

FIG. 6 is a diagram for explaining the manner in 
which the center of the expectation window is deter- 45 
mined by the first order line and field extrapolation 
methods. 

FIG. 7 is a diagram showing the memory configura- 
tion for storing the left-border coordinates of two fields. 

FIGS. 8A and 8B are a flow chart for calculating the 50 
center of the .expectation window by the line and field 
extrapolation methods. 

FIG. 9 is a diagram for explaining the manner in 
which three sample points are defined for the determi- 
nation of the reference level for the next line for the left 55 
border. 

FIG. 10 is a block diagram for explaining the manner 
in which the center of the expectation window is calcu- 
lated by the line and field extrapolation methods, as well 
as the sample and expectation points LSPL1 and LEPL, 60 
respectively for the left border. 

FIG. 11 is a block diagram for explaining the manner 
in which the positions of the sample and expectation 
points LSPL1 and LEPL, respectively are determined 
along a video scan line for the left border. 65 

FIG. 12 is a block diagram for explaining the manner 
in which the positions of the sample points LSPL2 and 
LSPR, and the end point LEPR of the expectation 
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window are determined along a video scan line for the 
left border. 

FIG. 13 is a diagram for explaining the manner in 
which the expectation and sample points for the right 
border are defined during the aortic valve simulation 
period. 

FIG. 14 illustrates the line periodic border model for 
the left border. 

FIG. 15 is a detailed circuit diagram of the dynamic 
reference level analog computer circuit. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

Referring now to the drawings, FIG. 1 illustrates, in 
block form, the real-time contour detector and data 
acquisition system for the left ventricular outline of the 
present invention. 

FIG. 2 illustrates the contour of the left ventricular 
outline as determined with the real-time contour detec- 
tor and data acquisition system of the invention. The 
X-ray image 11 of the contrast material filled left ventri- 
cle is stored for subsequent processing on film (cinean- 
giogram), video tape, or disc. The cineangiogram is 
converted into video format using a television camera 
13. The video signal from either the TV camera, video 
tape, or disc is the input signal to the system. The re- 
quired synchronization signals are generated by a sync 
generator 15. Before any contour algorithm is applied, 
several pre-processing operations are performed on the 
input signal to provide for a video signal with an im- 
proved signal to noise ratio and a restored DC level. 
The signal is amplified, and sync pulses which are pres- 
ent in the video signal are stripped off in an amplifier 
and sync stripper 17. Next, a 1 MHz, low-pass Bessel 
function filter 19 suppresses the upper frequencies 
which do not contain relevant contour information. 
The filtered signal is then clamped for DC level restora- 
tion in a video clamp 21. For display and visual feed- 
back, the input video signal is also passed through a 
delay line 23 which equalizes the total delay of the 
low-pass filter and the logic circuitry involved in the 
contour detection. The delayed signal is then applied to 
a summation amplifier 25 where the detected contour 
pulses are mixed into the video signal with the correct 
time relation. The total video signal is the output signal 
of the system and is used to display the X-ray image 
with the detected border on the TV monitor 27. 

The rest of the system is divided into two almost 
identical subsystems, one for the left border and one for 
the right border. Referring to FIG. 2, in this context, 
left and right are defined with respect to the chord from 
the left hand side of the designated aortic valve 29 plane 
to the apex 31 of the left ventricle, as viewed on the 
screen of the TV monitor 27. 

The left border subsystem 33 will now be described. 
An analog comparator 35 has a control input fed by a 
two-input OR gate 37. A starting point generator 39 is 
connected to one input of the OR gate and a digital 
processor unit 41 is connected to the other input of the 
OR gate. The analog comparator 35 is enabled for the 
first time during a video field at a starting point 43 
generated by the starting point generator 39. This start- 
ing point 43 can be manually positioned by the operator 
arbitrarily over the screen of the TV monitor 27, but is 
usually positioned at the aortic valve 29 plane. As soon 
as the video signal which is fed to the analog compara- 
tor 35 from the video clamp 21 reaches a reference level 
supplied from the dynamic reference level analog com- 
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puter circuit 45 and preset by the operator, the first left 
border point 47 is detected. As the comparator 35 
changes state, a narrow pulse is fed to the summation 
amplifier 25, where it is mixed into the video signal for 
visual feedback purposes. 5 

The position of a border point is preferably deter- 
mined in a 12-bit format by counting clock pulses from 
the rising edge of a horizontal sync pulse to this border 
point for the scan lines traversing the left ventricle; this 
is handled by the digital processor unit 41 which com- 10 
prises the circuitry depicted in FIGS. 10, 11 and 12. 
This digital processor unit 41 then determines the cen- 
ter, as well as the beginning — and end point of the new 
expectation window for the second left border point by 
an algorithm based on the location of the preceding 15 
border point and on the location of the corresponding 
border points in the preceding one or two video fields 
and also determines the location of three sample points 
used for the calculation of the reference level for the 
next border point. After the expectation window and 20 
sample points have been calculated, the left border 
RAM 49 is updated at the appropriate address with the 
last detected border point. To determine the location of 
sample points and expectation window along a video 
scan line the count is compared with the appropriate 25 
data words supplied by the digital processor. Three 
track and hold circuits 51 each sample and hold for a 
certain time interval on a line by line basis the voltage 
level at a respective one of the three sample points 
defined by the digital processor unit 41. The set of volt- 30 
age levels of the sample points is fed to the dynamic 
reference level analog computer circuit 45 which oper- 
ates on them to generate, in accordance with another 
algorithm, a new reference level for the comparator 35, 
resulting in the detection of the second left border point 35 
within the defined expectation window. For the third 
and following left border points, the digital processor 
unit 41 can apply each of the four possible combinations 
of line and field extrapolation principles for the calcula- 
tion of the new expectation window. 40 

The right border subsystem 53 will now be described. 
With a short fixed delay after the first left border point 
47 has been detected, an electronic generator 55 is 
started, simulating the aortic valve 29. Incorporation of 
the aortic valve simulation generator 55 is necessary 45 
because the automatic detection of this part of the ven- 
tricular outline is impossible due to the flow of contrast 
material into the aorta. The slope of the line represent- 
ing the aortic valve 29 is also manually adjustable. Vi- 
sual feedback resulting from the display of the X-ray 50 
picture on the screen of the TV monitor 27 together 
with the contour pulses, assists the operator in the oper- 
ation of the device. When the “artificial valve” meets 
the right border, its generation is terminated while the 
expectation window for the first actual right border 55 
point is determined. An analog comparator 57 has a 
control input fed by the digital processor unit 59. The 
aortic valve simulation generator 55 is connected to the 
digital processor unit 59. Digital processor unit 59 is 
identical to digital processor unit 41 except that unit 59 60 
has one additional input (on B-bus 85) to accommodate 
the output of aortic valve simulation generator 55. Dur- 
ing the aortic valve simulation period the digital proces- 
sor unit 59 determines the center of the expectation 
window from the position of the preceding border point 65 
in the current field plus a constant term which deter- 
mines the slope of the simulated aortic valve. The refer- 
ence level analog computer 61 applies a dc preset refer- 
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ence level to the analog comparator 57 during this per- 
iod. If the analog comparator 57 does not encounter a 
point with the same brightness level as the preset refer- 
ence level within the expectation window, then the end 
of the expectation window is defined as the contour 
point. However, if a point occurs with the correct refer- 
ence level within the defined expectation window, then 
this point is defined as the first right border point 63 and 
the aortic valve simulation generator 55 is stopped. As 
the comparator 57 changes state a narrow pulse is fed to 
the summation amplifier, where it is mixed into the 
video signal for visual feedback purposes. Because of 
the sudden change in direction from the simulated valve 
to the actual right border, the digital processor unit 59 
defines the center of the expectation window for the 
second right border point as having the same horizontal 
coordinate as the first right border point and also deter- 
mines the location of three sample points. After these 
coordinates have been calculated, the right border 
RAM 65 is updated at the appropriate address with the 
last detected border point. 

Three track and hold circuits 67 each sample and 
hold for a certain time interval on a line by line basis the 
voltage level at a respective one of the three sample 
points defined by the digital processor unit 59. The set 
of voltage levels of the sample points is fed to a dynamic 
reference level analog computer circuit 61, which oper- 
ates on them to generate, in accordance with another 
algorithm, a new reference level for the comparator 57. 
As soon as the video signal which is fed to the analog 
comparator 57 from the video clamp 21 reaches the 
reference level calculated by the dynamic reference 
level analog computer circuit 61, the second right bor- 
der point 69 is detected. For the third and following 
right border points, the digital processor unit 59 can 
apply each of the four possible combinations of line and 
field extrapolation principles for the calculation of the 
center of the expectation window. 

In the case where the video signal does not cross the 
reference level within the expectation window, forced 
border points are generated. This occurrence may be 
due to a sudden change in the video level at the border 
caused by intervening structures such as ribs or the 
diaphragm, or due to a sudden position change in the 
actual left ventricular border from previous border 
points. The possible situations are illustrated in FIG. 3 
for the left border. FIG. 3(a) shows the normal situation 
for a border point that is detected within the expecta- 
tion window. Consider now the case, as illustrated in 
FIG. 3 (b), in which the actual contour has shifted to the 
left such that the video signal is already above the refer- 
ence level at the beginning of the expectation window. 
A border point is then generated at the beginning of the 
expectation window as the comparator is enabled. FIG. 
3(c) shows the case in which the actual border has 
shifted to the right such that the video signal is below 
the reference level during the entire expectation win- 
dow period. This situation is recognized and a forced 
border point is then generated at the end of the expecta- 
tion window. 

The system is implemented so that no more than two 
contour points will occur on a video line traversing the 
left ventricle, one for the left border and one for the 
right border. 

The contour detection is stopped and an outpulse 71 
is generated by an outpulse circuit 73 connected to the 
computer interface 75 and the digital processor units 41, 
59 as soon as the distance between the left and right 
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border points is smaller than the adjustable ventricular 
apex 31 difference or if the left-hand side of the expecta- 
tion window for the right border occurs before a left 
border point has been detected, thus avoiding continua- 
tion of the detection process for the remaining portion 5 
of the field. 

The obtained border positions are transmitted from 
the analog comparators 35, 57 to a computer interface 
75 of conventional design. The computer interface de- 
termines the coordinates of the left ventricular border 10 
as detected by the contour detector in an 8-bit format 
over a number of video fields selected by the operator 
and controls their transfer to the core memory of a 
digital computer 77, thereby permitting real time opera- 
tion. For purposes of calculating size and shape of the 15 
left ventricle, it is sufficiently accurate to store each 
^-coordinate in an 8-bit format. The digital computer 
can be programmed to calculate the left ventricular 
volume from the detected border points using a known 
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If field m is an even field, h mI {(n+ 1)T] is defined as 

h m A(n+ DTI = ,)[(«+ 1)T] - g {m - 2) [{n+ 1)T] (4) 

The zeroth order field extrapolation principle takes 
only the last field into account. Applying both the line 
and field extrapolation principles simultaneously, the 
center of the expectation window ^ W [(«+1)T] on line 
(«+ 1) in the current field m is then defined as the aver- 
age position of the centers of the expectation windows 
determined separately by the line and field extrapola- 
tion methods. This results in: 

Field m is odd field: 

h m [(n + )7] = (5) 

2g m {nT) - g m [(n - 1)7] + 2g (m _ l) (nT) - g (m ^ 2) [(n + 1)7] 

2 

Field m is even field: 


h m [(n + 1)7] = 


2g m (nT) — g m [{n — 1)T] + 2g (m _ l) [(n + 1)7] — g( m _ 2 )[(n + 1)7] 
2 


( 6 ) 


mathematical model. 25 

The algorithm for determining the center of the ex- 
pectation window will now be discussed with reference 
to FIGS. 4, 5, and 6. A border point is denoted as 
g m (nT), where m refers to the field number and n to the 
line number. Using a first order line extrapolation prin- 30 
ciple, the center of the expectation window 
hmd( n + 1)T] on line («+l) is derived from the previ- 
ously detected border positions g m [( n ~ 1)T] and g m (nT) 
on video lines (n— 1) and n, respectively. The new posi- 
tion is linearly extrapolated, as shown in FIG. 4. The 35 
position of h mL [(n -f- 1)T] is calculated as 

h m J(*+ 1)T] = 2g m {nT) - g m [(n- 1)T] (1) 

The distance from the left hand side (LEPL) of the ^ 
expectation window to the center is denoted LEWL, 
the distance from the center to the right hand side 
(LEPR) is denoted LEWR, and T is the video scan line 
period. Right and left window widths are separately 
and manually adjustable. Using this linear extrapolation, 45 
the center of the window is generally taken as a good 
approximation of the next border point. 

Under low contrast conditions, it is sometimes advan- 
tageous not to extrapolate but rather to assign the center 
of the expectation window on line ( n + 1) the same x- 5Q 
coordinate as the last detected border point on line «, as 
illustrated in FIG. 5. This is referred to as the zeroth- 
order line extrapolation principle, and can be written as 

h mL [{n+l)T] = g>T) (2) 55 

Because the border position changes very little from 
field to field, the center of the expectation window in 
the current field can also be approximated by linearly 
extrapolating corresponding border points in the previ- 
ous two fields. This is referred to as the first order field 60 
extrapolation principle. 

To calculate the center of the expectation window 
using the first order field extrapolation principle, one 
must distinguish between even and odd fields. Assum- 
ing field m is an odd field, the center of the expectation 65 
window is defined as 

<wK«+l> n = 2g (m _i ) (nT) - g (m _ 2) [(n+l)T] (3) 


The zeroth order field extrapolation principle can be 
applied by substituting g (m _ 0 («T) and g (m _ {) [(«+ 1)T] 
in equations (5) and (6) respectively, instead of 
£(m-2)[(^ + 1)T]. 

The determination of the center of the expectation 
window using the first-order line and field extrapolation 
principles is illustrated in FIG. 6. The center of the 
expectation window can thus be determined in four 
different ways by selecting independently the zeroth- 
and first-order line and field extrapolation principles, 
respectively. 

A random access field memory 49 is included in the 
left border subsystem for the storage of the 12-bit coor- 
dinates, plus a control bit, from the last two fields which 
have been processed. An identical memory 65 is in- 
cluded in the right border subsystem. The memories are 
updated with each newly detected border point in the 
current field, and feed the associated digital processor 
unit 41 or 59 for determination of the center of the 
expectation window by the field extrapolation method. 
The memory configuration for storing the left border 
coordinates is illustrated in FIG. 7. With each coordi- 
nate available in a 12-bit format, the thirteenth bit is 
used as a control bit, which is zero, if there is no border 
point on a scan line and one if there is a border point. 
Storing the 13 bits present at the DATA IN input at a 
location specified by a 10-bit address whose eight least 
significant bits are determined by counting the scan 
lines during a video field, whose 9th bit is always zero, 
and whose most significant bit can be 0 or 1 depending 
on whether one wants to access the coordinates of an 
odd or even field respectively, is accomplished by hav- 
ing the READ/WRITE and CONTROL ENABLE 
inputs both low. Reading data from a specified memory 
location occurs when the READ/WRITE input is high 
and the CONTROL ENABLE low. After the DATA 
OUT outputs have stabilized, these data can be loaded 
into the memory buffer register 79. Using tri-state buff- 
ers 81, this information can be applied to the associated 
digital processor unit during the appropriate time peri- 
ods to calculate the center of the expectation window. 

FIG. 8 is a flow chart of the procedure used to calcu- 
late the center of the expectation window for the left 
border on line («+ 1) in field m. The processing cycle is 
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divided into 11 time states by digital processor unit 41 
and the time states are designated tsO through tslO. 

The field extrapolation principle is enabled by acti- 
vating a switch shown at the top of the flow diagram. 

So long as this switch has not been activated, the center 5 
of the expectation window is determined using the line 
extrapolation principle only (this will be clarified later 
on). To always start in the same phase after having 
activated the switch, the first complete odd field after 
this moment is defined as the first field with FE (field io 
enable) = 1. The most significant address bit A9 is made 
low at the beginning of each video line in an odd field to 
access the lower half of the memory 49 and is made high 
at the beginning of each video line in an even field. The 
first full line in a field after the occurrence of the verti- 15 
cal synchronization pulse is assigned the value n = L 
The calculation of the center of the expectation win- 
dow /i m [(«+l)T] is only carried out if a left-border 
point has been detected at the current line n . If no left- 
border point has been detected during a scan line, then 20 
g m (nT) = 0 is stored into memory 49 with a correspond- 
ing control bit equal to zero, indicating that all zeros 
have been stored at this location. As a result, all scan 
lines not traversing the left ventricle are assigned a 
control bit equal to zero. Because of the vertical move- 25 
ment of the entire left ventricle over a heart cycle plus 
its contraction and ejection movements, it is important 
to know if actual contour points were detected at the 
selected video lines in the previous fields. 

Assuming a left-border point has been detected, the 30 
most significant address bit is inverted at the beginning 
of time-state tsO to access the border coordinate in the 
previous field. Also, a 1 is added to the address if the 
current field is an even field (compare eqs. (3) and (4)). 
During time state tsl, the center of the expectation 35 
window is calculated according to the line extrapola- 
tion principle as 2gJjiT)-g m [(n — 1)T]. This calculation 
is done in the digital processor 41. The result is stored in 
a 16 bit register 87 (FIG. 10) of the same type as used for 
the memory buffer register 79. This 16-bit register is 40 
referred to as ACL (accumulator left border) and its 
contents as M(ACL); therefore, 

M(ACL) = 2g>T)-U(«-l)T] 

45 

If the zeroth-order line extrapolation is selected, g m (riT) 
is substituted for g m [(n— 1)T] and the result is 

M(ACL) - IgJnT) - gJnT) = g m (nT) 
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end of ts3, which yields M(MBR): = 2M(MBR). These 
new data are added to the contents of the accumulator 
during ts4. Assuming that one had selected the first- 
order line extrapolation principle and activated the field 
extrapolation switch, the contents of ACL would be 

M(ACL) - 2£ m (nT>g m {(rt — 1)T] +2g {m „ 1)(Nl T) 

At the beginning of ts3, the address is set for the reading 
of the data from field (m— 2) during a later time state. 
The most significant address bit A9 is again inverted, 
thus giving it the initial value it had at the beginning of 
the video line, and a 1 is added to the current line num- 
ber. The next step during ts5 is to shift the data in the 
MBR one position to the left, resulting again in 

M(MBR) — £( m _i } (N 1 T) 

Meanwhile, the coordinate g( m - 2 )[( w + l)T] is read 
from memory 49, but not yet loaded into the MBR. At 
the end of ts5, the system checks whether the zeroth- 
order field extrapolation principle was selected and 
whether the control bit for g( m „ 2) [(rc-i-l)T] is equal to 
zero. If the answer to one or both of these checks is true, 
then the zeroth-order field extrapolation will be ap- 
plied. The border coordinate g (m _ 2) [(« + 1)T] will not be 
stored in the MBR, but the old data g (m _ 1} ( N{T) are 
subtracted from the accumulator contents during ts7. If 
this is not the case, then g( W - 2 )[(« + 1)T] is stored in the 
MBR at the end of ts2 and the contents of MBR are 
subtracted from the accumulator contents. 

Earlier, the conditions were described for which the 
MBR was cleared for the remainder of the processing 
cycle (FF = 0). It is clear that, under these conditions, 
the contents of MBR do not change even if one wants to 
load g( m _ 2 )[(« + l)T], because MBR remains cleared 
independently of the clock pulse. 

These same conditions are also important for the next 
step. If FF=0, the accumulator contains the center of 
the expectation window calculated according to the 
line extrapolation principle only. For this condition, the 
clock to the accumulator is inhibited during ts8 and 
nothing happens. If FF ^ 0, then the accumulator con- 
tains the sum of the centers of the expectation windows 
as calculated according to the line and field extrapola- 
tion principles separately. To obtain the average of 
these centers, the contents of the accumulator are 
shifted one position left during ts8, resulting in 


The field memory 49 is enabled at the beginning of 
time state tsl and, with the correct address applied, 
&(m-i)(NiT) is read and stored into the memory buffer 
register at the end of ts2, that is 

M(MBR) - g {m _ X) { N,T) 

where Nj = n if the current field is an odd field and Nj 
= n + 1 for an even field. At the end of ts2, it is also 
checked if the field extrapolation principle was not 
selected (FE=0?) and if the control bit for g^-i^NjT) 60 
is equal to zero. If the answer to one or both of these 
questions is true, then the MBR is cleared for the re- 
mainder of this processing cycle. This situation is re- 
ferred to as FF:=0. This means that the center of the 
expectation window will be determined only by the line 65 
extrapolation principle. Independent of whether the 
data in the MBR are unequal to zero or not, the contents 
of MBR are shifted over one position to the right at the 


M(ACL) - h m [(n+\)T) 

At the beginning of ts7, the extra 1 in the memory ad- 
dress is removed and the last detected border point 
g m (nT) is stored in memory during ts8 and ts9. 

The expectation window width LEWL is subtracted 
from A m [(«+1)T] during ts9 and the result is stored in 
the accumulator. These data are also stored in D flip 
flop registers; the contents of these registers then repre- 
sent the left-hand side LEPL of the expectation win- 
dow. Similarly, the sample window width LSWL is 
subtracted from the contents of the accumulator during 
tslO and stored in ACL at the end of tslO. The final 
result is then that the accumulator contains the position 
of LSPL1 for the next video line in the current field. 

The processing cycle is now complete and we return 
to the beginning of the flow chart for the next process- 
ing cycle. If line n was not the last line in the video field, 
then the address is incremented by one and bit A9 is 


55 
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reset to zero if the current field is an odd field and set to 
one for an even field. With n being the last line in the 
field, the line counter is reset during the vertical syn- 
chronization pulse and the whole cycle is repeated. 

A digital processor unit 41, 59 is included in the sys- 5 
tem for the calculation of the center, the beginning and 
end points of the expectation window and the positions 
of the sample points. 

Formulas (5) and (6) are given for the general case 
that the center of the expectation window is determined 10 
using both the first-order line and field extrapolation 
principle. The zeroth-order line extrapolation principle 
can be applied by substituting g m {riT) for g m [(« — 1)T] at 
the time of calculation. Similarly, the zeroth-order field 
extrapolation principle can be applied by substituting I 5 
g{m- i)(«T) and g (m _ i)[(«+l)T] in equations (5) and (6) 
respectively, for g( m « 2 )[(»+ 1)T]. With each border 
position known in a 12-bit format, the center of the 
expectation window is also determined in a 12-bit for- 
mat, applying formulas (5) and (6). However, we also 20 
must determine the beginning and end points of the 
expectation window to enable the comparator, and the 
sample points to calculate the dynamic reference level. 

If it is assumed for simplicity, that only the first-order 
line extrapolation principle is applied, the beginning and 25 
end points of the expectation window on line ( n + 1) are 
shown in FIG. 4 for the left border. These points are 
denoted LEPL and LEPR, respectively. According to 
the way the system has been implemented, these begin- 
ning and end points are defined as 30 

LEPL = A m [(w+ 1)T] - LEWL (7) 


and 


LEPR = LEPL + LEWL + LEWR 


35 

( 8 ) 


that is, LEPR is determined from the position of LEPL. 

Three sample points are defined for the left border as 
illustrated in FIG. 9. The sample points are denoted 40 
LSPR, LSPL1 and LSPL2, respectively. Assuming the 
last detected border point is gJjiT) on line n , the video 
signal is sampled on the same line inside the ventricle at 
a distance LSWR from g m (nT) and on the next line 
(h + 1) in the background area at distances LSWL and 45 
LSWL/2 before the left-hand side of the expectation 
window LEPL, respectively. The video samples at the 
three defined sample points are used to calculate the 
reference level for line («-f 1)- FIG. 9 clearly shows that 
the position of the sample points LSPL1 and LSPL2 are 50 
determined with respect to the center of the expectation 
window, thereby taking into account the local direction 
and movement of the left ventricular outline. This re- 
sults in a more accurate determination of the reference 
level for the next border point. 55 

In formulas: 

Sample point LSPR is determined from the detected 
border point g m (nT) on line n as 

LSPR = g m (nT) + LSWR (9) ^ 

and sample point LSPL1 is calculated from LEPL as 
LSPLl = LEPL - LSWL (10) 

In the implemented system, sample point LSPL2 is 65 
derived from LSPLl as 

(ID 


The selected window widths LEWL, LEWR, LSWL 
and LSWR are available in binary format. 

The digital processor unit 41 comprises the subsys- 
tems shown in FIGS. 10, 11 and 12 and functions to 
determine the center, the beginning and end points of 
the expectation window and the positions of the sample 
points. These subsystems are described in more detail in 
the following paragraphs. 

The digital processor unit 41 contains an arithmetic 
unit as shown schematically in FIG. 10 for the left bor- 
der. 

Using tri-state buffers, the addends are all applied to 
the A-bus 83, and the augends or subtrahends to the 
B-bus 85. The result of an arithmetic addition or sub- 
traction is stored in the accumulator ACL 87 and can be 
applied again to the A-bus for the next arithmetic opera- 
tion. 

After the left border point g m {nT) has been detected, 
its 12-bit coordinate is available at the inputs of the 
register 89 and buffer 91. Consider the case when the 
zeroth-order line extrapolation principle is selected. At 
the start of time-state tsl, the 12 bits are stored in the 
register 89. With the tri-state buffers 91 and 93 enabled 
during tsl, 2 g m (nT) is applied to the A-bus 83 and 
g m (n T) to the B-bus 85 (the multiplication by two has 
been achieved by shifting the buffer output one bit posi- 
tion on the bus). The Arithmetic Logic Unit 95 (ALU) 
control being set for subtraction, the resulting 2 gJnT) 
- T) is stored at the end of tsl into the accumulator 
ACL 87. 

For the first-order line extrapolation principle, 2 
g m {n T) is again applied to the A-bus 83 during tsl 
through buffer 91. However, at this moment, g m {nT) is 
not stored in the register 89; during tsl the information 
in this register is then still g m [(« — 1)T] from the previ- 
ous line and this coordinate is applied to the B-bus 85. 
At the end of tsl, the resulting 2 gJnT) - g m [(n — 1)T] is 
stored in the accumulator 87. 

The contents of the accumulator 87 can be applied to 
the A-bus 83 at the appropriate time states as defined in 
FIG. 8 for additional additions and subtractions. The 
data from the left border RAM 49 is applied to the 
B-bus through the buffer 97. 

The binary word representing LEWL is applied to 
buffer 99 and the binary word representing LSWL to 
buffer 101. This allows the calculation of LSPLl and 
LEPL. After the position of LEPL has been calculated 
it is stored in register 103. At the end of the processing 
cycle, the position of LSPLl is stored in the accumula- 
tor 87. 

The position determination of these points on a scan 
line is also performed in the digital processor unit and is 
described with the block diagram of FIG. 11. The data 
stored in the ACL 87 and register 103, representing 
LSPLl and LEPL respectively, are connected to the 
multiplexer 105. From the beginning of a scan line until 
the position of LSPLl has been determined, the multi- 
plexer control is set such that the 12 bits representing 
LSPLl are applied to one side of the 12 bit comparator 
107. The other side of the comparator is connected to 
the outputs of a 50 MHz counter 109. This counter 109 
starts counting the clock pulses at the end of the hori- 
zontal sync pulse for the scan lines traversing the left 
ventricle. As soon as the contents of the counter 109 
equal the 12-bit data word for LSPLl; the output of the 
comparator 107 becomes logically high, indicating the 
position of LSPLl. The multiplexer control is now 


LSPL2 = LSPLl + LSWL/2 
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changed such that the 12 bits for LEPL are applied to 
the digital comparator 107, whose output becomes logi- 
cally low again. The counter 109 continues to count and 
the output of the comparator 107 becomes high again as 
the state of the counter 109 equals the 12-bit data word 5 
for LEPL. This initiates the enabling of the analog 
comparator 35. As the video signal crosses the refer- 
ence level, the analog comparator 35 changes state and 
the counter 109 stops counting clock pulses. The state 
of the counter 109 thus represents the new border coor- 10 
dinate on this particular line. 

The position determination of the points LSPL2, 
LEPR and LSPR are also performed in the digital pro- 
cessor unit and are described with the block diagrams in 
FIG. 12. For the position determination of LSPL2, the 15 
4 bits representing LS WL/2 are applied to one side of a 
4 bit digital comparator 111 and the outputs from a 4-bit 
counter 113 to the other side. This 50 MHz counter is 
enabled as soon as LSPL1 has been detected. The 
counter 113 then starts counting until its state equals 
LSWL/2, so that the output of the comparator 111 
becommes high indicating the position of LSPL2. The 
position determination of the points LEPR and LSPR 
are performed in identical ways. For LEPR, a 5-bit ^ 
counter 115, connected to a 5-bit comparator 117, is 2 
enabled by LEPL and the counter 115 keeps on count- 
ing until its state equals LEWL + LEWR, which is 
applied to the other side of the comparator 117. If the 
analog comparator 35 did not detect a border point ^ 
within the expectation window, then the comparator 35 
is disabled at this moment, resulting in the generation of 
a forced border point at the end of the expectation 
window. For LSPR, a 4-bit counter 119, connected to a 
4-bit comparator 121, is enabled by g m (nT) — that is the 35 
detected border point on this particular line — and the 
counter keeps on counting until its state equals LSWR, 
which is applied to the other side of the comparator 121. 

The way the expectation and sample points for the 
right border are determined in the digital processor unit ^ 
59 are equivalent to the way described above for the left 
border. However, for the right border we also have to 
include the aortic valve generation, so that the proces- 
sor unit includes one more input of B-bus 85 (FIG. 10) 
to accommodate the output of aortic valve simulation 45 
generator 55. The equations for the right border points 
are described below. 

The positions of the expectation and sample points 
during the aortic valve simulation are given in FIG. 13. 

So long as this simulated valve plane is generated, the 50 
field extrapolation principle is inhibited and only the 
zeroth-order line extrapolation principle is applied, nec- 
essary for the generation of a straight line. 

The beginning point REPL of the expectation win- 
dow on line («+ 1) is determined by 55 

REPL = gJnT) tAV- REWL (12) 

where Av is the aortic valve width, given in an 8-bit 
format. The end point of the expectation window dur- ^ 
ing the aortic valve plane generation is defined as 

REPR = REPL + REWL (13) 
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Sample points RSPR1 and RSPR2 on line n are derived 
from gJjiT) as 

RSPRl - g m {nV) + RSWR/2 (15) 

and 

RSPR2 = gJnV) + RSWR (16) 

The general formulas for the expectation and sample 
points after the aortic valve simulation period can be 
derived from formulas (12), (13) and (14). For the gen- 
eral case with the center of the expectation window 
denoted h m [(n + 1)T], equation (12) becomes 

REPL = «J(*+1)T] + AV - REWL (17) 

The end point REPR is defined as 

REPR = REPL + REWL + REWR (18) 

Equation (14) can now be written as 

RSPL « A ffl ((n+1)T] + AV - REWL - RSWL (19) 

The formulas for RSPRl and RSPR2 remain, of course, 
unchanged. In order that the timing for the right border 
arithmetic unit remains unchanged, the term AV is 
always included, but is made equal to zero outside the 
aortic valve simulation period. 

Next the algorithm for dynamically adjusting the 
reference level according to local brightness levels on a 
line by line basis will be discussed. The brightness level 
changes along a border because of shading, nonhomo- 
geneous distribution of the contrast agent in the left 
ventricle, and overlapping of Roentgen shadows from 
other organs and structures, such as diaphragms and 
ribs. Three sample points have been defined in FIG. 9 
for the left border. 

The reference level for the border point on line 
(n + 1) is then calculated as 

, \ (20) 
V K ffl + 1) = f J+ V c 

where a is a proportionality factor, V c is a constant 
voltage level, VR is the video sample at LSPR, and 
VL1 and VL2 the video samples at LSPL1 and LSPL2, 
respectively. By taking the average value of VL1 and 
VL2, either side of the contour has the same weight. 
The function of the two sample points at the back- 
ground area will be explained later. 

The behavior of the detection process can be ana- 
lyzed by assuming a line periodic video signal, as shown 
in FIG. 14 for the left border. In order to simplify the 
calculations, one left sample poinit LSPL is assumed and 
both sample points LSPL and LSPR have the same 
distance A to the border point. With the slope of the 
video signal within the ventricle denoted by y l V/s and 
in the background area by y 2 V/s (with y 2 <7i), and, 
assuming b/y x < BP(n) < b/y v 4- A, the final value of 
the reference level can be calculated as 


Sample point RSPL is defined as 

RSPL - gJnX) + AV - REWL - RSWL = 
REPL - RSWL 


65 
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( 21 ) 


( 14 ) 
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A border point is said to be acceptable if it occurs on 
the ventricular slope within a distance A from the knee 
in the border model* This results in one sample point in 
the background area and the other on the ventricular 
slope. 5 

The requirement for the reference level is then 

b < V w/ < b + A* (22) 

Solving (21) and (22) yields « 0 


Y fri - y *> + b 


< a < 1 - 


K 

b + Ay! 


(23) 


A practical value for a is 0.8 to 0.9, and for the distance 1 
A, approximately 0.5 us. 

Using the linearized border model of FIG. 14 it can 
be theoretically derived that the optimum relation be- 
tween a and V c under varying parameters y u y 2 , and b ^ 

1S v c =yi -a) (24) 


where b MIN is the lowest knee level in the image. Ap- 
proximating b MIN in the current field with the measured 
(VL1/2 -J- VL2/2) M1N of the previous field and apply- 25 
ing equation (24) then only requires an adjustment of 
the a factor for the appropriate V re f> as determined 
according to equation (21). 

The two sample points LSPL1 and LSPL2 are impor- 
tant in areas of low contrast where the actual video 30 
signal may be different from the assumed border model. 

At those places, it often occurs that the slope y 2 be- 
comes negative, resulting in a too high reference level, 
as calculated from equation (20). This situation is char- 
acterized by VL1>VL2 and is recognized in the sys- 35 
tern. Under this condition, an error term V erwr = VL1/2 
- VL2/2 is determined and subtracted from the term 
VL1/2 + VL2/2 in the formula for the reference level. 
The effect is that the reference level is determined as if 
the border signal had a background slope of y 2 = 0. On 40 
these bases, the final formula for the reference level can 
be written as 




a_( ILL VL2 
2 \ 2 ^ 2 


- Kr„r+ VR 



(25) 

45 


where 


VIA 


VL2 


, if 


VL\ 


VL2 


= 0 , 


.„ VL 1 „ VL2 
lf _<_ 
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Referring to FIG. 15, there is shown a detailed circuit 55 
diagram of the dynamic reference level analog com- 
puter circuit 45 for operating on the sample points 
LSPR, LSP1, and LSPL2 to generate a reference level 
for the analog comparator 35. 

The video samples V R , VL1/2, and VL2/2, and the 60 
negative error signal are applied to the negative input of 
the operational amplifier 123 which is used as a —1 
amplifier. The signals VL1/2 and VL2/2 at half the 
amplitude of the original signal can be obtained by use 
of an accurate voltage divider (not shown). The output 65 
of amplifier 123 is the negative of the sum of these sig- 
nals. The error term is or is not included, depending on 
the setting of the switch 125. The a/2 factor in the 
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formula for the reference level is implemented with a 
resistive divider 127. The voltage follower operational 
amplifier 129 acts as a buffer and its output signal equals 
the product of a/2 and the output of amplifier 123 
where 0 < a < 1. The negative adjustable level of 
— V c volts is obtained with a potentiometer 131 con- 
nected to a zener diode 133 stabilized reference voltage. 
The — V c volts and the output signal from amplifier 129 
are applied to the negative input of the operational 
amplifier 135 which is used as a —1 amplifier. The 
output of this amplifier then equals the negative of the 
sum of the latter two signals which is the required dy- 
namic reference level. This signal is connected to one of 
the inputs of a single-pole, double-throw junction FET 
switch 137. Depending on the state of the switch, the 
output is then either the dynamic reference level or the 
output signal from the operational amplifier 139, which 
provides an adjustable DC level. This DC level is ap- 
plied when a constant reference level is selected for the 
border detection and for the detection of the first left- 
border point when the dynamic reference level mode is 
selected. The desired switch setting is achieved by ap- 
plying the appropriate logic signal at the input of the 
high-speed switch driver 141. 

The error term is generated with the operational 
amplifiers 143 and 145. The sampled video level VL1/2 
is applied to the positive input of the operational ampli- 
fier 143 and VL2/2, to the negative input. The opera- 
tional amplifier 143 is used as a difference amplifier, 
resulting in an output voltage equal to the difference of 
the latter two signals. When (VL1/2) exceeds (VL2/2) 
is determined by use of the precision limiter circuit 147. 
With the positive input of amplifier 145 connected to 
ground, the output of the precision limiter is just the 
required error signal. 

Obviously, numerous modifications and variations of 
the present invention are possible in light of the above 
teachings. It is therefore to be understood that within 
the scope of the appended claims, the invention may be 
practiced otherwise than as specifically described 
herein. 

What is claimed as new and desired to be secured by 
letters patent of the United States is: 

1. A real-time contour detector and data acquisition 
system for an angiographic apparatus having video scan 
means for converting X-ray images of a cardiac cham- 
ber characterized by a change in brightness level com- 
pared with its surrounding into video format and dis- 
playing the X-ray images in recurring video fields and 
comprising: 

sampling means adapted to be coupled to the video 
scan means for sampling the brightness level of 
predetermined sample points of the current video 
field in the region of the last detected border point; 

reference level means connected to the sampling 
means and responsive to the sampling of the bright- 
ness level of the predetermined sample points for 
generating a reference brightness level at which 
the next border point can be detected; 

comparator means connected to the reference level 
means for comparing the brightness level along a 
video scan line in the current video field with that 
of the reference brightness level to detect the next 
border point; 

field memory means for storing the coordinates of 
detected border points from two preceding video 
fields; 
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processor means connected to the comparator means, 
the field memory means, and the sampling means 
for enabling the comparator means only within a 
narrow expectation window on the next video scan 
line centered about the expected next border posi- 5 
tion and for determining the sample points to be 
sampled by the sampling means; and 
data acquisition means connected to the comparator 
means for acquiring the detected border points. 

2. The system recited in claim 1 wherein the refer- 10 
ence level means generates a reference level representa- 
tive of 

f( i r L + i l i - >*)+>', 

where 

a is a proportionality factor, 

V c is a constant voltage level, 

VR is a voltage sample corresponding to a point 
within said cardiac chamber, 

VL1 and VL2 are voltage samples corresponding to 
points in said surrounding, and 

v KL1 VL2 . r VL\ VL2 

V error — 2 2 * U 2 > 2 

= o, 

3. The system described in claim 1 wherein the refer- 
ence level means generates a reference level which is a 
function of the brightness of a sample point within said 
cardiac chamber and the brightness at two sample 


18 

points outside said cardiac chamber, the position of said 
outside sample points being established with respect to 
the center of said expectation window so as to take into 
account the local direction and movement of the car- 
diac chamber contour. 

4. The system as set forth in claim 1 wherein means 
are included for generating a forced border point if the 
brightness level encountered during the period of the 
expectation window fails to exceed said reference 
brightness level. 

5. The system recited in claim 1 wherein means are 
included for generating a forced border point when the 
brightness level encountered at the onset of the expecta- 
tion window exceeds said reference brightness level. 

6. A system as claimed in claim 1 wherein said proces- 
sor means determines the coordinates of the narrow 
expectation window for the next border point as a func- 
tion of the location of the detected border points of at 

2 q least one preceding video line in the current video field. 

7. A system as recited in claim 1 wherein said proces- 
sor means determines the coordinates of the narrow 
expectation window for the next border point as a func- 
tion of the location of the corresponding border point in 

25 at least one preceding video field and as a function of 
the location of the detected border point in at least one 
preceding video line of the current video field. 

8. The system set forth in claim 1 wherein said pro- 
cessor means determines the coordinates of a set of 

30 sample points as a function of the location of the last 
detected border point and the location of the center of 
the narrow expectation window for the next border 
point. 
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